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We present quantum calculations of the relaxation matrix for the Q branch of N 2 at room tempera- 
ture using a recently proposed N 2 -N 2 rigid rotor potential. Close coupling calculations were comple- 
mented by coupled states studies at high energies and provide about 10200 two-body state-to state 
cross sections from which the needed one-body cross-sections may be obtained. For such temper- 
atures, convergence has to be thoroughly analyzed since such conditions are close to the limit of 
current computational feasibility. This has been done using complementary calculations based on the 
energy corrected sudden formalism. Agreement of these quantum predictions with experimental data 
is good, but the main goal of this work is to provide a benchmark relaxation matrix for testing more 
approximate methods which remain of a great utility for complex molecular systems at room (and 
higher) temperatures. © 2014 AIP Publishing LLC. [http://dx.doi.Org/10.1063/l.4862082] 


I. INTRODUCTION 

Collision-induced rotational energy transfer in molecular 
gases plays an important role in a variety of gas phase pro- 
cesses, such as pressure broadening and shifting and/or mix- 
ing of spectral lines.' The most accurate calculation of the 
corresponding cross-sections involves solving sets of coupled 
differential equations using the close coupling (CC) method.” 
However, even with present day computers, such calculations 
remain prohibitively expensive for complex molecular sys- 
tems (e.g., a mixture of polyatomic molecules). As a result 
of these limitations, various approximate models have been 
developed. They include fitting laws and more accurate scal- 
ing approaches; amongst the latter, those derived from the 
energy corrected sudden approximation (ECS) seem to be 
very powerful.' Alternatively, some researchers have devel- 
oped purely classical^ or semi-classical"'’^ calculations of the 
corresponding relaxation matrices. 

In a recent paper,® some of us have shown that it is possi- 
ble to introduce line coupling effects within the semi-classical 
Robert-Bonamy formalism,^ leading to a new method for cal- 
culating the off-diagonal elements of the relaxation matrix. 
We have also shown in that paper that considering line cou- 
pling leads to better agreement between semi-classical and 
fully quantum linewidths (diagonal elements of the relaxation 
matrix W), starting from the same intermolecular potential for 
the N 2 -N 2 pair. That system was chosen for benchmarking 
purposes: Thibault et had previously reported detailed 
fully quantum calculations of the half-widths of the isotropic 
Raman Q lines, based on a very recent potential energy sur- 
face (PES) proposed by Gomez et al.^^ Before investigating 
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the validity of this new semi-classical approach for the cal- 
culation of the off-diagonal elements of W, it is necessary to 
dispose of fully quantum non-diagonal elements which will 
be used as a reference for testing the semi-classical results. 
Such CC calculations were previously reported by Green and 
Huo," but based on the PES of van der Avoird,'^ which has 
been shown to need some refinements.^’®’'^ More recently, 
Eonfria et al . have reported CC two-body cross-sections at 
total energies too low to obtain converged results for temper- 
atures higher than 50 K. It is therefore necessary to consider- 
ably extend these calculations in order to work at room tem- 
perature. This is the purpose of the present paper; the compar- 
ison with the semi-classical results obtained from the method 
of Ref. 6 will be presented in a forthcoming paper. 

In Sec. II, an outline of the quantum formalism and nu- 
merical calculations is given together with a brief summary of 
the ECS formalism. Section III compares the current results 
with experimental data and previous calculations and provides 
a brief discussion of the convergence of the calculations, be- 
fore a brief summary of our findings. 

II. THEORY 

A. Some definitions 

Within the impact and binary collisions approximations, 
the effects of molecular collisions on isotropic Raman spec- 
tra is described by a relaxation matrix Wij'\, j\, T) whose 
rows and columns are labeled by spectral lines'’'®’'® (for 
simplicity, 7 1 means here Q(/i); j\ and f\ refer to the spec- 
trally active molecule). At low pressures, only diagonal 
elements are important since they define the half-widths 
of non-overlapping lines (y( 7 , T)=W ( 7 , 7 ; T)). At higher 
pressures, off-diagonal elements become important and 
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describe line mixing transfer of intensity between overlap- 
ping lines. When the PES does not contain any vibrational 
dependence, W is real and may be written in terms of one- 
body cross-sections which are, for off-diagonal elements, the 
negative of ordinary state-to-state cross-sections for rotational 
excitation (or deexcitation): 




ll2V 

2nc 


cr\j[JuT) 


fl2V 

2nc 


o-(7i 




( 1 ) 

where «2 is the perturber density and v is the mean relative 
velocity. These cross-sections may be written as a weighted 
sum of two body rotational state to state cross-sections: 


o-(yi ^ y[; t’) = X! X! ^ j'l’ t), (2) 

h A 


where ]2 and i '2 refer to the bath molecule and p(/2) is the nor- 
malized rotational population of level 72 (including the spin 
factor). Finally, we recall that the two-body cross-sections, 
at a given temperature T, are obtained as a thermal average 
over the initial relative kinetic energy of inelastic cross- 
sections cr(/i,j2 ^/l,/ 2; EkinY 

o-(;i , h j[ , i 2 \ T) = 7— —2 / o’(;'i , h j[ . 72; Eki„) 

{K-bI ) Je, 

X exp( Ekin / ^BE)Ekin^ Eklfi , (3) 


where £s is the minimum kinetic energy for the levels / 1 and 
j '2 to become accessible. 

In cases where the collision dynamics do not depend on 
vibrational motion, a sum rule exists which says that, for each 
line, the diagonal elements of W are equal to the negative of 
the sum of the off-diagonal elements:'® 

(4a) 


Y. oih ^ j[\ T) = a" (;i, h-T) . (4b) 

Another important relation is the detailed balance rela- 
tionship: 

P O'l) o-(yi ^ j [ ; T) = p(;i')o-(;i' -> j \ ; T), (5) 

which results, at a more elementary level, from the detailed 
balance for the two-body cross-sections: 

P0'l)p0'2)cr(yi, 72 ^ j[, j' 2 , T) 

= p(7i)p(7>Op;2^;i,72;E), (6) 

the latter resulting from the microscopic reversibility, at a 
more elementary level: 

(2Ji + 1) (272 -b 1) EkO- (71, 72 ^ j'l, j' 2 , Ek) 

^{2j[+\){2j'^+\){Ek+E{j,)-E{j[)+E{j2)-E{j'^)) 

X j '2 7i> 72; Ek+E{jx)-E{j[)+E{j 2 )-E{j 2 )). 

(7) 


B. Quantum dynamical calculations 

Quantum dynamical calculations were performed on the 
four dimensional ab initio potential energy surface of Ref. 10 
determined by using symmetry adapted perturbation theory. 
This interaction potential, between rigid rotors, is expected 
to be the most accurate to date.^"*"’*^ In particular, it pro- 
vides accurate self-broadening coefficients over a wide range 
of temperatures (77-1700 K) for Raman isotropic Q lines.*’® 

Close-coupling and coupled states calculations were per- 
formed with the MOLSCAT code and its parallel version.** 
Technical details handling the PES and the dynamical calcu- 
lations can be found in Refs. 8 and 9. Both variants of N2 
have been considered, ortho (0N2; even j's) and para (pN2; 
odd fs). The PES does not allow interconversion between 
ortho and para species. Following previous work,"’*"*’*® the 
calculations were performed assuming that the two collid- 
ing molecules are distinguishable, even for 0N2-0N2 or pN2- 
pN2 collisions, thus neglecting (small) resonance-exchange 
effects. Moreover, since the experimental results we are com- 
paring with*® involve vibrational excitation of one of the two 
N2 molecules, such a contribution will be much smaller than 
in the case of two molecules in the same vibrational state. Dy- 
namical calculations were performed at various total energies 
Ej, and thus at various kinetic energies Eki„ — Ej — E^otiji, 
72) for a given total internal rotational energy given by E^otiji, 
72) = Bjiiji 4- 1) -b 572(72 + 1) with the rotational constant 
B — 1.998 cm“*. Here we only present results for an ortho 
active molecule ij\,j'\ even). 

Two-body state-to-state cross-sections o{j\j 2 fij' 2 ', 
Ekin) for ortho-para nitrogen collisions were calculated us- 
ing the CC method over a grid of 150 total energies up to 450 
cm“*, thus including the (asymptotically) open two-body ro- 
tational energy levels lower than Erot(€>, 13) = 447.6 cm“*, 
and completed by 15 coupled states calculations for total en- 
ergies between 500 and 1500 cm“* (thus including the open 
two-body rotational levels up to 71 = 22,72 = 15). Therefore 
two-body rate constants (see below) involving rotational en- 
ergy levels lower than £^((6, 13) are CC/CS (coupled states) 
results while results for rotational energy levels higher are 
fully CS. 

Similarly, calculations for 0N2-0N2 were performed us- 
ing the CC method over a grid of 254 total energies up to 
505 cm“*, including the open two-body rotational energy lev- 
els lower than £ror(6, 14) = 503.55 cm“*, and completed by 
15 coupled states calculations up to 1500 cm“* (including 
two-body rotational levels up to 71 = 24, 72 = 12). There- 
fore two-body rate constants for rotational energy levels lower 
than Ei-ot(6, 14) are CC/CS results while results for rotational 
energy levels higher are fully CS. In order to illustrate the 
progress of the computational power, note that the CC cal- 
culations of Ref. 1 1 were limited to total energies up to 200 
cm“* while CS calculations were obtained up to 700 cm“*. 

Because of the limitations of our grid in kinetic 
energies, for 0N2-PN2 and 0N2-0N2 we have only re- 
tained the cross-sections (j(jij 2 —>■ /1/2; Tj, respec- 
tively, between the very first 70 two-body rotational 
levels (i.e., up to 71 = 10,72 = 15 associated with a rotational 
energy of about 700 cm“*) and those between the very first 
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75 two-body rotational levels (i.e., up to j\ = 2,72 = 18 asso- 
ciated with a rotational energy of about 695 cm“'). However 
experiments measure effective one particle cross sections'"^’ 
a{ji f\\ T) that are obtained from Eq. (2). In evaluating 
this double sum, it is necessary to include all the signihcantly 
populated bath level 72 , as well as all the significant two-body 
cross sections (i.e., those with a small inelasticity). The most 
populated level at 300 K corresponds to 7 ss 6 - 8 . We include 
all the levels up to 72 =18, which corresponds to 98% of the 
population and the sum over i '2 incorporates contributions up 
to i '2 — ji + 20. Therefore, due to the limits of our CC/CS 
calculations, it was necessary to check the convergence of our 
CC/CS one body cross-sections. We have therefore developed 
a method inspired by the previous work of Green and Huo:^' 
we have used the ECS formalism to generate, when necessary, 
extrapolated two-body cross-sections (i.e., those not obtained 
from our CC/CS calculations). If the ECS contribution to 
Eq. (2) remains small, then the one body cross-section may 
be reasonably considered as a CC/CS result; otherwise it will 
not be retained as a reference value. It is therefore necessary 
to recall briefly some basic ECS relations allowing the calcu- 
lation of all the two body state to state cross sections from a 
limited set of “fundamental” ones: 


Q(Ly,L2) = o{LuL2^0,Q\T). ( 8 ) 

Note here that our definition of these basic cross-sections 
is different from that of Green and Huo"’ who used Q^^(L \ , 
L 2 ) = cr( 0 , 0 — ^ Li, L 2 ), so that the two sets are related by the 
detailed balance Eq. ( 6 ). 


C. ECS formalism 

Within the approach of Green,'® a non-diagonal two- 
body cross-section can be deduced from the following scaling 
relation: 


h ^ j[, 72 ; T) = ( 27 ; + 1)(272 + 1) 


^(2Li + 1)(2L2+1) 

Li,L2 



Li 

0 



L2 

0 



X 


24 + (f2 (LuL2)tcf 
24-t- (f2(7i, j2)rcf 


Q(LuL2), 


(9) 


where ( 1 1 1) is a 3j-symbol. Eq. (9) is applied only for ener- 
getically downward collisions (positive transfer from rotation 
to translation) while upward cross-sections are obtained from 
detailed balance (Eq. (6)). According to De Pristo etal.,^^ the 
adiabaticity factor [ ^24+^(j'i’h)T^ which accounts for the 
non-resonant character of the collisions is defined in terms 
of an average collision duration — U/v where is an ad- 
justable scaling length. The frequency factor represents the 
dominant inelasticity for level 7172 and it can be defined in 
various ways. Eollowing the previous analysis of Green,'® we 
used: 


( 10 ) 


III. RESULTS AND DISCUSSION 

In the following, results will be indifferently presented as 
relaxation matrix elements Wj'j in units of 10 “^ cm“* atm~' 
(Eq. (1 )) or as one body excitation rates ri 2 Va{j\ — >• 7 '[ ; T), in 
units of /xs“' Torr”* which were adopted by Sitz and Farrow*® 
to present their experimental data. 

A. Determination of the ECS scaling length 

As may be seen from Eq. ( 8 ), the basic Q{L \ , L 2 ) are also 
given by the CC/CS calculations since they simply correspond 
to the usual inelastic cross-sections. They are available up to 
(Li = 18; Z .2 = 18) and are given in Table I. Therefore the next 
step is the determination of the only unknown parameter: 
This was determined via a least-square fit of Eq. (9) to the 
complete set of 10232 available CC/CS cross-sections. The 
best fit was obtained for Z^. = 2.4 A. Among the 10 232 cross 
sections, 72% of the ECS predictions were in error by less 
than a factor of 2 , with a root mean square relative error of 
0 . 12 ; 28% were in error by more, with an extremely large rms 
relative error of 2.3, reflecting very large errors in a very small 
number of tiny cross-sections, mostly corresponding to very 
large inelasticity. 

Before continuing, it is worthwhile comparing these re- 
sults with those previously obtained by Huo and Green (ESQ 
method in Ref. 1 la)) with a slightly different procedure : they 
fitted both the scaling length and the fundamental , L 2 ) 

rates, while we keep these fixed to the CC/CS values. They 
obtained a scaling length of 2.5 A, very close to our value. 
Now, by comparing the last column of Table I, to the basic 
fLse(0, 0 — >• Li, L 2 ) determined by Huo and Green (column 
labeled LSQ in Table III of Ref. 11a)) it appears that their 
fLse(0, 0 — >• Li, L 2 ) are close to our CC/CS values (deduced 
from the Q{L\, L 2 ), via the detailed balance), at least for the 
lowest {L\,L 2 ) couples. It is therefore not surprising that the 
two methods converge towards very similar solutions for the 
effective one-body rate constants. 

B. Determination of the <t(/i y'i ; T) cross-sections 

On this basis it was possible to calculate one body cross- 
sections in different ways, from Eq. (2). Some results are 
given in Table II. Results labeled CC/CS used only available 
CC/CS two body cross-sections, setting all missing ones to 
zero. Results labeled ECS correspond to the same sum except 
that we used ECS two body cross-sections. Comparison of 
these two columns gives an idea of the accuracy of the ECS 
predictions for the (highly averaged) one body cross-sections. 
The third column gives an ECS estimation of the missing 
CC/CS terms in the double sum of Eq. (2) (those set to zero). 
Finally, the CC/CS results were supplemented by adding the 
ECS correction to the CC/CS results, leading to the final re- 
sults (last column). 

When the ECS correction is much smaller than the pure 
CC/CS contribution (less than about a few per cent), the fi- 
nal result may be considered as a confident result: “quite a 
pure CC/CS cross-section.” This is the case for the cr(/j =4 
^ /i; T — 298 K) cross-sections appearing in Table 11(a). 


^ (yj> ji) — 2 fi( 7 i -t- 72 )- 
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TABLE I. Basic rates for N2-N2 treated as distinguishable particles at 
298 K, in /xs“^ Torr“^. Column (1) corresponds to the a{L\, L 2 — ^ 00) as 
given by our CC/CS calculations. Column (2) gives the corresponding cr(00 
^ Li, L 2 ) obtained from column (1) via the detailed balance (Eq. (6)). To be 
compared with column labeled LSQ in Table III of Ref. 1 1 (a). 


Li 

L2 

(1) 

(2) 

0 

2 

4.26 X 10“' 

2.011 

0 

4 

6.14 X 10^2 

0.456 

0 

6 

2.4 X 10^2 

0.208 

0 

8 

1.182 X 10^2 

0.1 

0 

10 

6.39 X 10-3 

4.67 X 10-3 

0 

12 

3.98 X 10-3 

2.22 X 10-3 

0 

14 

2.42 X 10-3 

9.35 X 10-3 

0 

16 

1.32 X 10-3 

3.2 X 10-3 

0 

18 

6.79 X 10-^ 

9.4 X 10-^ 

2 

2 

1.65 X 10-* 

3.673 

2 

4 

3.651 X 10-3 

1.28 

2 

6 

1.11 X 10-3 

0.456 

2 

8 

5.4 X 10-3 

0.217 

2 

10 

2.72 X 10-3 

9.38 X 10-3 

2 

12 

1.55 X 10-3 

4.09 X 10-3 

2 

14 

8.49 X 10-^ 

1.55 X 10-3 

2 

16 

4.46 X 10-^ 

5.09 X 10-3 

2 

18 

2.24 X lO-'* 

1.46 X 10-3 

4 

4 

2.324 X 10-3 

1.282 

4 

6 

8.51 X 10-3 

5.49 X 10-* 

4 

8 

3.52 X 10-3 

2.23 X 10-* 

4 

10 

2.09 X 10-3 

1.13 X 10-* 

4 

12 

1.074 X 10-3 

4.454 X 10-3 

4 

14 

5.72 X 10-^ 

1.64 X 10-3 

4 

16 

3.02 X 10-^ 

5.42 X 10-3 

4 

18 

1.58 X lO-'* 

1.62 X 10-3 

6 

6 

8.14 X 10-3 

6.14 X 10-* 

6 

8 

3.84 X 10-3 

2.84 X 10-* 

6 

10 

1.45 X 10-3 

9.22 X 10-3 

6 

12 

7.11 X 10-^ 

3.45 X 10-3 

6 

14 

3.59 X 10-^ 

1.2 X 10-3 

6 

16 

1.94 X 10-^ 

4.07 X 10-3 

6 

18 

9.64 X 10-3 

1.16 X 10-3 

8 

8 

3.464 X 10-3 

2.51 X 10-* 

8 

10 

1.185 X 10-3 

7.365 X 10-3 

8 

12 

3.7 X lO-'* 

1.76 X 10-3 

8 

14 

2.29 X 10-^ 

7.52 X 10-3 

8 

16 

1.154 X 10-^ 

2.38 X 10-3 

8 

18 

5.82 X 10-3 

6.86 X 10-^ 

10 

10 

9.19 X 10-^ 

4.9 X 10-3 

10 

12 

2.76 X lO-'* 

1.125 X 10-3 

10 

14 

1.04 X lO-'* 

2.92 X 10-3 

10 

16 

5.29 X 10-3 

9.34 X 10-^ 

10 

18 

3.57 X 10-3 

3.61 X 10-"* 

12 

12 

2.76 X lO-'* 

8.61 X 10-3 

12 

14 

7.4 X 10-3 

1.6 X 10-3 

12 

16 

4.48 X 10-3 

6.05 X 10-^ 

12 

18 

2.14 X 10-3 

1.66 X 10-^ 

14 

14 

2.125 X 10-3 

3.16 X lO-'* 

14 

16 

1.86 X 10-3 

1.74 X 10-^ 

14 

18 

1.3 X 10-3 

6.93 X 10-3 

16 

16 

1.71 X 10-3 

1.0 X 10-^ 

16 

18 

7.9 X 10-3 

2.64 X 10-3 

18 

18 

4.8 X 10-3 

9.21 X 10-3 


TABLE II. One body rate constants for 71 ^ /i transitions at 298 K, in 
10“^ cm~^ atm~^, calculated from Eq. (2). Column labeled CC/CS: sum of 
the available CC/CS two body rates. Column labeled Equivalent ECS: sum 
of the same two body rates but obtained from the ECS scaling relation (Eq. 
(9)). Column labeled ECS complement: sum of the two body rates missing 
in the CC/CS calculation and calculated from the ECS scaling relation. Final 
result corresponds to column (2) + column (4). 


/i 

CC/CS 

Equivalent ECS ECS complement 

Final result 

0 

2.02 

(a) For;i = 4 

1.88 

0.04 

2.06 

2 

11.43 

10.34 

0.18 

11.61 

6 

14.85 

12.0 

0.24 

15.09 

8 

8.37 

6.17 

0.18 

8.55 

10 

4.5 

3.3 

0.15 

4.65 

12 

2.12 

1.54 

0.12 

2.24 

14 

0.8 

0.58 

0.083 

0.883 

16 

0.2 

0.17 

0.05 

0.25 

18 

0.008 

0.045 

0.037 

0.045 

0 

0.168 

(b) For ji = 14 
0.098 

0.016 

0.184 

2 

0.873 

0.55 

0.09 

0.963 

4 

1.55 

1.11 

0.18 

1.73 

6 

2.3 

1.79 

0.26 

2.56 

8 

3.21 

2.54 

0.35 

3.56 

10 

5.0 

3.34 

0.47 

5.47 

12 

10.6 

4.93 

0.74 

11.34 

16 

3.13 

2.78 

0.81 


18 

0.1 

0.99 

0.79 



First note the very good accuracy of the ECS predictions for 
small inelasticities. Consider for instance the case ji = 4 
^ /i = 6 and the sum (ji = 4, j 2 j[ = 6, T) 

h 

plotted as function of 72 in Fig. 1(a). As it appears the most 
important differences between CC/CS and ECS results appear 
mainly for high 72 values which are weighted by small popu- 
lations in Eq. (2). Consequently, as shown in Fig. 2(a), the two 
calculations of Eq. (2) lead to values differing only by 20%. 
One may also observe the good convergence of the CC/CS 
sum. 

Things are more complex for higher 71 values, as may be 
seen from Table 11(b), which corresponds to 71 = 14. As ex- 
pected, the accuracy of the ECS predictions is not as good for 
such a high 71 value. Moreover, while the ECS supplement 
remains small (when compared to the CC/CS sum) up to /i 
= 12 , it becomes important for f\ — 16 — 18 (and higher /i) 
and is probably underestimated. This may be easily under- 
stood from a similar analysis as above. We consider the case 
7 i = 14 — >• /i = 16. As shown in Fig. 1(b) important differ- 
ences exist for the sum Y ^ (/i = 14, 72 ^ j[ — 16, j'2, T) 
J2 

calculated by CC/CS or ECS methods. The convergence of 
the summation (Eq. (2)) over 72 values is shown in Fig. 2(b). 
The CC/CS and ECS calculations start to diverge for 72 > 8 
but that divergence stops abruptly for 72 > 12 due to the lack 
of CC/CS values. Therefore an important ECS complement is 
needed in that case as it appears in Table 11(b). Hence a(ji 
= 14 ^ /i > 14; r = 298 K) cannot be considered as 
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FIG. 1. (,/l. J 2 — ^ j[, 72’ ^ function of j 2 '. 

22 

(a) fory'i = 4 ;/i = 6 and (b) for^'i = 14 ;/, = 16. 

benchmark CC/CS values. As a consequence, the linewidth 
cross section cr*^(14, 14; T = 298 K) as deduced from the sum 
rule (Eq. (4b)) cannot be labeled as a CC/CS result. 

C. Comparison with experimentai data 

As mentioned previously, rotational excitation rates for 
N 2 (v = 1)-N2(v = 0) collisions have been measured by Sitz 
and Farrow'® at 298 K by a pump-probe technique and may 
be used to check the consistency of our final results. Indeed 
as shown by Huo and Green," rotational collisional rates 
do not strongly depend on the different vibrational states of 
the two colliding partners. As can be seen from Table III, 
there is a good agreement overall between experiment and 
theory: most of our CC/CS final results are inside or very 
close to the error bars of the experimental data. Unfortu- 
nately, there are large uncertainties on some of the experi- 
mental results and further experimental investigation would 
be of great interest before any new analysis of the remain- 
ing discrepancy between theory and experiment. As outlined 
previously,"’'"' some of the differences may be due to inac- 
curacy in the PES and/or to the neglect of exchange contribu- 


72 max 

FIG. 2. Convergence of the sum cr(yi ^ 7j;7) = ^ piji) 

72=0 

as function of Amax: (a) for ji = 4, /i =6 
■ii 

and (b) for ji = 14, /i = 16. 

tions. While interesting, such problems are beyond the scope 
of the present paper. Table IV summarizes the main result of 
the present work: the fully quantum relaxation matrix at room 
temperature for the isotropic Raman Q branch of N 2 . Note 
that the diagonal elements are those of Ref. 9. Indeed, since 
the rule in the present work is to only consider “quite pure 
CC/CS results” it is not possible to obtain fully converged 
linewidths from Eq. (4b) for all Q(ji) lines. The convergence 
of the diagonal elements was considered in a different way in 
Ref. 9, by summing first over the post-collisional rotational 
quantum numbers (j'l, j' 2 ), analyzing at that step the con- 
vergence and performing then the thermal average and fi- 
nally summing over j 2 the partial pressure broadening cross- 
sections a(ji,j 2 ', T) weighted by the populations p(/ 2 )- Never- 
theless, the collisional halfwidths obtained summing the off- 
diagonal elements of Table IV (Eq. (4b)) reproduce the main 
trend with y'l of the diagonal elements® and lead to very close 
values to the latter at least for the first few values. Due to 
the lack of off-diagonal terms for jiifi) >12 the agreement 
strongly deteriorates for ji > 10. 
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TABLE III. One body rate constants for ji f i transitions at 298 K, in 

/ns“* Torr^*. Comparison of our CC/CS results, completed by ECS correc- 
tions, with the experimental data of Sitz and Farrow.*^ 


h 

/i 

Experiment 

CC/CS final results 

0 

2 

6.64 ± 1.18 

7.23 

0 

4 

3.76 ± 0.83 

3.78 

0 

6 

2.73 ± 0.61 

2.58 

0 

8 

0.86 ±0.51 

1.69 

0 

10 

0.64 ±0.12 

0.94 

0 

12 

0.29 ± 0.06 

0.45 

0 

14 

0.22 ± 0.04 

0.175 

2 

4 

5.13 ±0.58 

4.52 

2 

6 

2.40 ± 0.44 

2.83 

2 

8 

1.52 ±0.34 

1.81 

2 

10 

0.97 ± 0.22 

1.01 

2 

12 

0.28 ±0.14 

0.49 

2 

14 

0.12 ±0.04 

0.19 

4 

6 

4.7 ± 0.6 

3.74 

4 

8 

2.2 ± 0.4 

2.12 

4 

10 

1.4 ±0.2 

1.15 

4 

12 

0.54 ±0.13 

0.56 

4 

14 

0.21 ± 0.04 

0.22 

6 

8 

3.37 ± 0.47 

3.18 

6 

10 

2.22 ± 0.29 

1.51 

6 

12 

0.71 ±0.14 

0.77 

6 

14 

0.23 ± 0.05 

0.28 

8 

10 

2.52 ± 0.41 

2.65 

8 

12 

1.12 ±0.22 

1.07 

8 

14 

0.29 ± 0.07 

0.39 

10 

12 

2.68 ± 0.43 

2.26 

10 

14 

1.04 ±0.13 

0.72 

12 

14 

1.83 ±0.26 

1.96 


TABLE IV. CC/CS relaxation matrix IV for N 2 -N 2 at 298 K (The Wyj 
elements are expressed in 10“^ cm“* atm“*). Diagonal elements have been 
reported in Ref. 9. 


j 

r 

0 

2 

4 

6 

8 

10 

12 

14 

0 

68.2 

-6.2 

-2.1 

- 1.22 

-0.82 

-0.54 

-0.34 

-0.185 

2 - 

-29.2 

50.3 - 

- 11.6 

-6.3 

-4.1 

-2.7 

- 1.69 

- 0.966 

4 - 

-15.3 - 

- 18.3 

46.0 

- 13.0 

-7.53 

-4.8 

-3.05 

- 1.745 

6 - 

- 10.46 - 

-11.5 - 

- 15.1 

44.4 - 

- 13.14 

-7.3 

-4.9 

-2.575 

8 

-6.84 

-7.33 

-8.6 

- 12.9 

44.0 - 

- 12.57 

-6.67 

-3.59 

10 

-3.81 

-4.1 

-4.67 

-6.12 - 

- 10.73 

42.3 - 

-11.8 

-5.5 

12 

- 1.84 

- 1.97 

-2.26 

-3.13 

-4.34 

-9.13 

40.1 

-11.37 

14 

-0.7 

-0.77 

-0.88 

- 1.13 

- 1.6 

-2.9 

-7.92 

34.1 


IV. CONCLUSION 

To summarize, using the most recent rigid rotor N 2 -N 2 
potential reported so far, we have carried out, at the CC and 
CS levels, room temperature calculations of the relaxation 
matrix for the isotropic Raman Q branch of N 2 . Such a tem- 
perature gives rise numerical problems associated with the 
convergence of the quantum calculations. However, we have 
shown that it is possible to quantify the convergence by us- 
ing the ECS approximation. We then demonstrated that for 
moderate values of the rotational quantum number j < 12, 
the convergence of our CC/CS cross sections has been prac- 
tically achieved, since ECS complements are limited to a few 
percent. This fully quantum relaxation matrix is expected to 
be a benchmark for testing the validity of more approximate 
methods using the same PES, in particular the semi-classical 
approach proposed in Ref. 6. 
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